Zero Energy Solutions and Vortices in Schrodinger Equations 
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Two-dimensional Schrodinger equations with rotationaUy symmetric potentials (Va(p) — 
— a^5ra/9^'^'"^^'with p = ^y x'^ -\- y'^ and a 7^ 0) are shown to have zero energy states contained in 
conjugate spaces of Gel'fand triplets. For the zero energy eigenvalue the equations for all a are 
reduced to the same equation representing two-dimensional free motions in the constant potential 
Va = —Qa in terms of the conformal mappings of C,a ~ with z = x -\- iy. Namely, the zero energy 
eigenstates are described by the plane waves with the fixed wave numbers ka ~ y/mga/{h/2n) in 
the mapped spaces. All the zero energy states are infinitely degenerate as same as the case of the 
parabolic potential barrier (PPB) shown in ref. 8. Following hydrodynamical arguments, we see 
that such states describe stationary flows round the origin, which are represented by the complex 
velocity potentials W = Paz'^ , ( Pa being a complex number) and their linear combinations create 
almost arbitrary vortex patterns. Examples of the vortex patterns in constant potntials and PPB 
are presented. We see that any vortices cannot be created in the linear combinations of the plane- 
wave solutions and the degenerate states which are not described by plane- waves play essential roles 
to create vortices. In the extension to three-dimensional problems with potentials being separable 
into 2 + 1 dimensions we show that the states in three dimensions have the same structure as the 
two-dimensional states with the zero energy but they can generally have non-zero total energies. 
PACS numbers: 03.65.-w, 03.50.-z 



I. INTRODUCTION 

It is known that scattering states and unstable states 
like resonances are generally described by states in con- 
jugate spaces of Gel'fand triplets (GT) An example 
of such states for the parabolic potential barrier (PPB) 
V = — m7^a;^/2 in one dimension has been studied by 
many authors [2-7]. It has been shown that the one- 
dimensional PPB has pure imaginary energy eigenvalues 
=Fi(n-|-l/2)?i7 with n = 0, 1, 2, • • • , and the eigenfunctions 
are generalized functions in the conjugate space 5(M)^ 
of GT described by 5(R) C ^^(M) c 5(M) ^ , where 5(R) 
and £^(R), respectively, stand for a Schwartz space and 
a Lebesgue space In general the energy eigenval- 

ues £ of the conjugate spaces in GT are expressed by 
pairs of complex conjugates such that £ = e =F 17 with 
e,7 G M and the states with the =F sign, respectively, 
represent resonanc-decay and resonance-formation pro- 
cesses. This pairing property of the energy eigenvalues 
indicates that states in higher-dimensional PPB possi- 
bly have zero energy eigenstates. In fact the PPB in 
two dimensions (generally in even dimensions) has zero 
energy eigenvalue, which is included in the eigenvalues 
expressed by ^iirix — ny)h'y. We see that the zero en- 
ergy states are obtained for zero and positive integers 
satisfying Ux = riy and then they are infinitely degener- 
ate. The zero energy states are interpreted as stationary 
flows round the centre of the PPB Furthermore, fol- 
lowing hydrodynamics, it has been also shown that some 
of such flows can be expressed by complex velocity poten- 
tials and various vortex structures appear in the linear 
combinations of the infinitely degenerate states. Con- 
sidering that states in the conjugate spaces of GT are 
generally not normalizable but currents of those states 



are observable in quantum mechanics, the quantities ob- 
served in physical processes should be based on the prob- 
ability currents such as currents in hydrodynamics. Hy- 
drodynamical considerations will play a very important 
role in the investigations of quantum physics in GT. Hy- 
drodynamical approach of quantum mechanics was vig- 
orously investigated in the earlier stage of the develop- 
ment of quantum mechanics [9-16]. Vortices were ex- 
tensively examined by Hirschfelder [17-20] and a review 
article was written by Ghosh and Deb ]^. It should 
be noted that such hydrodynamical idea is still useful in 
present-day quantum physics [p[ p^j2^ . Actually prob- 
lems of vortices appear in many aspect of present-day 
physics such as vortex matters (vortex lattices) p^ , p5[ , 
vortices in non-ncutral plasma [26-29] and Bose-Einstcin 
gases [30-34] and so on. The vortex problems will hold a 
very important position in the hydrodynamical approach 
of quantum mechanics. As noted in ref. j^, the station- 
ary flows in the two-dimensional PPB can create almost 
arbitrary patterns of vortices because of the infinite de- 
generacy. PPB potentials can be a good approximation 
to the repulsive forces that are very weak at the centre of 
the forces such that harmonic oscillators are a good ap- 
proximation to the attractive forces being very weak at 
the centre. In fact PPB has been applied in some chemi- 
cal problems [35-37] The PPB, however, is a very special 
potential and then it seems to be difficult that the results 
of PPB extend to more general potentials. 

In this article we shall investigate stationary flows in 
more general types of potentials from the hydrodynam- 
ical point of view. Especially, wc study stationary flows 
in two dimensions discussed in the PPB because in- 
teresting quantities in hydrodynamics such as complex 
velocity potentials and vortices are definable in the di- 
mensions. Furthermore it is well known that in two- 
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dimensional hydrodynamical problems conformal map- 
pings can be a very strong tool to investigate velocities, 
complex velocity potentials and vortices [38-41]. Some 
solutions solved in a special aspect are possibly extend 
to others in terms of conformal mappings. Particularly 
we will pay attention to the stationary flows that are de- 
scribed by the zero energy eigenstates in the PPB 
Such zero energy states can also play an interesting role 
in statistical mechanics in GT, where a new type of en- 
tropy arises from the freedom with respect to the imagi- 
nary parts of energy eigenvalues [42-44]. That is to say, 
even if a many-particle system is in the ground state with 
a fixed energy and then it has no freedom arising from 
the real energy eigenvalue, it can still have the freedom 
with respect to the imaginary parts. Remembering the 
fact that all the energy eigenvalues in GT are expressed 
in the pairs of complex conjugates, we can understand 
the situation very easily, because the imaginary part of 
every complex energy can be cancelled out by adding the 
conjugate energy in many-particle systems. An example 
was presented in ref. [Q for the PPB, where the burst 
of entropy from the new entropy was studied in thermal 
non-equilibrium. If we can find such zero energy states 
in more general types of potentials, we can investigate 
many-particle systems from a very new aspect where ev- 
ery state with a fixed energy can still have huge variety 
arising from the degeneracy of zero energy states. This 
article has two themes. One of the themes is to show 
the fact that rotational symmetric potentials of the type 
Vaip) = -a^5a/5^^°"^^/^ with p = ^/jF+y^ and a ^ 
have the same zero energy solutions as those obtained 
in the PPB in two dimensions [||. The other is vortex 
formation in terms of the zero energy solutions. We shall 
show that, as far as the zero energy solutions are con- 
cerned, Schrodinger equations with the rotational sym- 
metric potentials Va{p) for all values of a excepta = 
are reduced to the same equation by using conformal 
mappings. We shall also see that the conformal map- 
pings which are known to be very powerful tools in two- 
dimensional hydrodynamics become powerful tools also 
in the hydrodynamical approach of quantum mechanics 
and vortex patterns for all the potentials Va{p) can be 
investigated in a very simple method. 

We shall perform our considerations as following; in 
section 2 general property of Schrodinger equations with 
rotational symmetric potentials is investigated in terms 
of conformal mappings. In section 3 it is shown that, for 
zero energy solutions, all the equations in the mapped 
spaces can be reduced to the same equation describing 
free motions in a constant potential. This means that, 
as far as the zero energy eigenstates are concerned, all 
the symmetric potentials have the same solutions with 
the infinite degeneracy as those obtained in the PPB [^ . 
Boundary conditions are also discussed there. Follow- 
ing the considerations of the PPB, hydrodynamical argu- 
ments are performed and velocity, complex velocity po- 



tentials and vortices are investigated in section 4. In sec- 
tion 5 an extension of the argument to three-dimensional 
problems are carried out and vortices are studied in three 
dimensions. Non plane- wave solutions for the zero energy 
are briefly discussed in section 6. Remarks on non-zero 
energy solutions are presented in section 7. Some re- 
marks and comments on the present work are given in 
section 8. 

II. CONFORMAL MAPPINGS OF 
SCHRODINGER EQUATIONS WITH 
SYMMETRIC POTENTIALS 

We shall investigate the general structure of 
Schrodinger equations 

d 

where the Hamiltonian H is described by rotational sym- 
metric potentials in two-dimensional space {x,y). The 
eigenvalue problems with the energy eigenvalue £ are ex- 
plicitly written as 

[- — A + Va{p)]^{x,y)=£^j{x,y), (1) 
2m 

where a e M (a 7^ 0), 

dx^ dy"^ ' 

with p — x"^ + y^, m and ga are, respectively, the mass 
of the particle and the coupling constant. Note here that 
Va represents repulsive potentials for {ga > 0, a > 1) 
and {ga < 0, a < 1) and attractive potentials for {ga > 
0, a < 1) and {ga < 0, a > 1). Since we investigate 
the equations in the conjugate spaces of GT, the energy 
eigenvalues £ of (Q) are generally complex numbers. 

Following the hydrodynamical argument [38-41], let us 
consider the following conformal mappings; 

(^a — z°-, with z = x + iy. (2) 

Note that the conformal mappings are singular at the 
origin except in the cases for a ^positive integers, and the 
conformal mapping for a = 1 is trivial because nothing 
is changed by the mapping. We further notice that a 
complex factor A can be multiplied in the mappings such 
as C = Az"", which will be discussed in the case for A — 
e'" with a e K. When we use the notation 

we see that 

Ua = p"^ cos aip, Va — p"" sill aip, (3) 



2 



where Ua, Va G M. and p = \/ + y^, V' = arctan(2//x). 
Using the notations 

Pa = + and = av3, 

we have 

Ua = Pa COS ((fa and = Pa sin (4) 

In the (Ma,Wa) plane the equations (|l|) are written down 
as 

^2^2(a-l)/a[_^^^ _ VK, «a) = ^ V^K, I'a), (5) 

2 m 

where 

We can rewrite the equations as 
2 m 

(6) 

Exchanging the second term of the left-hand side and the 
term of the right-hand side, we obtain 

2 m 

(7) 

It is quite interesting that we can read this equations as 
follows; the eigenvalue problem for the potential Va{p) 
in the (x, y) plane given by (^) is replaced by the eigen- 
value problem for the potential Vifa{pa) in the (ua,fa) 
plane, where the roles of the eigenvalue £ and the cou- 
pling constant are exchanged. We may consider that 
this relation represents a kind of duality between the en- 
ergy and the coupling constant. From the relation we 
see that by solving the eigenvalues for the fixed £ in the 
(ua, Va) plane we can determine the strength of the cou- 
pling constant ga to reproduce the eigenvalue £ for the 
potential Va{p) in the {x, y) plane. We shall return to the 
relations between the problems for Va in the {x^y) plane 
and Vi/a{pa) in the {ua,Va) plane in section 7, because 
this theme is not the main subject of this section. 

Here let us briefly comment on the conformal mappings 
C^a = z"" ■ We see that the transformation maps the part 
of the (x, y) plane described by < p < oo, < < 7r/a 
on the upper half-plane of the {uajVa) plane for a > 
and the lower half-plane for a < 0. Note here that the 
maps on the part of the (ua,fa) plane with the angle 
LPa = — a can be carried out by using the conformal 
mappings 

Ca(«) - ^"e-^". (8) 



In the maps the variables 

Ua{a) = p° cos(a(^ — a) and Va{Q) = p° sm{ap — a) 

(9) 

should be used. We also have the relations 

Ua{a) — Ua cos a -\- V a sva. a and Wa(a) = fa cos a — Ua sin a. 

(10) 

Of course, the relations Wa(0) = Ua and fa(0) = Va are 
obvious. 

It will be better to comment on the meaning of the 
choice of the variables Ua and Va given in (3). It is obvi- 
ous that Ua and Va are not suitable variables to represent 
the states having definite properties with respect to rota- 
tions, such as the states with definite angular momentum, 
in comparison with the polar coordinates p and ip. In 
the following discussions, however, we will be interested 
only in the states describing stationary flows that are ba- 
sic elements in hydrodynamics. In general the stationary 
flows, such as those in scattering problems, cannot be de- 
scribed by the states with definite angular momentum, 
because every stationary flow has the specific directions 
representing the incoming and out-going flows. (Exam- 
ples of the stationary flows in PPB will be presented in 
section 4. See figs. 1 and 2.) Such stationary flows, 
of course, have no definite rotational symmetry except 
rotations with respect to some specific angles. We can 
understand such situations by considering the fact that 
the directions of the incoming flows are chosen by hand 
in scattering experiments. Actually it will be shown that 
the freedom of the phase a in the conformal mapping (8) 
is related to such choices. (See section 3.) The choice 
of the variables Ua and Va is, therefore, important in the 
following hydrodynamical approach, where the relations 
between the potentials Va{p) with different values of a 
are studied. An explicit example of the difference be- 
tween the choice of the polar coordinates and that of Ua 
and Va has been shown in the case of PPB in section 3 
ofref. [8]. 

III. ZERO ENERGY SOLUTIONS OF THE 
SCHRODINGER EQUATIONS 

We shall here study the special solutions having zero 
energy eigenvalue £ — Q. As noted in section 1, en- 
ergy eigenvalues in GT are generally complex and all en- 
ergy eigenvalues appear as pairs of complex conjugates 
like £ =F *7 (S: 7 G Hi) fill- This indicates that, provided 
that a potential in one dimensional space has pure imag- 
inary eigenvalues, the potentials extended in two dimen- 
sions possibly have zero energy states. This situation 
really occurs in parabolic potential barriers (PPB), that 
is, one dimensional PPB has pure imaginary eigenvalues 
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[2-7] and hence two dimensional (generally in even di- 
mensions) PPB has zero energy states that are described 
by the stationary flows round the origin and infinitely 
degenerate |^. Let us investigate zero energy solutions 
for the potentials Va{p)- 

A. Zero energy solutions 

We see that for the zero energy S = the Schrodinger 
equations obtained by the conformal mappings be- 
come very simple such that 

h^Aa -5a] V'K,Wa) = 0. (11) 

Note that the zero energy solutions have no time depen- 
dence. It is remarkable that the equation becomes same 
for all a, that is, the potential is expressed by the con- 
stant Qa for all a. As far as the zero energy solutions are 
concerned, the equations transformed by the conformal 
mappings can be written in the same form for all the 
potentials Va{p) with a ^ 0. It should be noticed here 
that only in the case of the constant potential Vi = —gi 
for a = 1 the energy eigenvalues can take arbitrary values 
satisfying th condition £+gi > 0, because the right-hand 
side of (6) does not have any p dependence. 

It is trivial that the solutions of (11) are given by the 
two-dimensional plane waves with the energy ga- The 
solutions are, therefore, represented by 

^±(pJ = 7V,e±*'="W•''^ for5, >0 (12) 

and 

</.±(pJ = Af,e±'="W''», forg„<0 (13) 

where the angle denotes the moving direction of the 
plane wave in the (ua, Va) plane, 

ka{0) = \/'2m\ga\/fi ■ (cos 0, sin 6*) and — {ua,Va) are 
two-dimensional vectors, and Na and Ma are in general 
complex numbers. Comparing the equation ka{0) ■ pa = 
ka {ua cos 9 ^ Va sin 6) where ka = ^/2m\ga\/ H with Ua (a) 
of (10), we see that the angle 9 can be adjusted to the 
phase a introduced in the conformal mappings (8). By 
using the phase a and the variable Ua{a) the solutions 
(12) and (13) are written by 

^o±(Ma(a)) ^iV.e^''^"""^"), for<?„>0 (14) 

and 

</)±(ua(a)) ^Mae*'^"""^"), iorga<0. (15) 

We shall use the representations given in (14) and (15) in 
the following discussions. Note here that, taking acount 
of the relations 

Uo(±7r/2) = iWa, Ua(±vr) -Uo, (16) 



we can represent all the solutions of (14) and (15) by 

tpQ {ua{a)) and (l)Q{ua{a)) with — tt < a < tt. 

(17) 

We also notice that the solutions il>^{ua{a)) for ga > 
are expressed by the plane waves with the fixed momen- 
tum Pa — \J2mga^ whereas those (/)^(ita(a)) for g^, < 
are expressed by exponential growing or dumping func- 
tions. This difference is essential, because the plane- 
wave solutions can always be the states contained in the 
conjugate spaces of GT of which nuclear space is given 
by Schwartz space whereas the exponential growing 
functions such as ef cos(ai^-a) ^\\\^ q <- cos(a</7 — a) can- 
not find a simple nuclear space for arbitrary values of a. 
From now on we shall mainly discuss on the plane-wave 
solutions for > 0. 

Let us summarize the main results in the cases of 
ga > 0. 

(1) All the potentials written by Va{p) have zero energy 
eigenstates in GT. 

(2) All the solutions with the zero energy can be ex- 
pressed by the plane-wave with the fixed momentum 
Pa = ^/2mga in the iua,Va) plane. 

(3) The zero energy solutions have an infinite freedom 
arising from the arbitrary angle — tt < a < tt that corre- 
sponds to the freedom of the angle between the incoming 
particle and the x axis, which is given by (tt — a) /a. 

(4) In the case of the constant potential corresponding to 
a = 1 , though we have the same solutions obtained in the 
above arguments, their energy eigenvalues need not equal 
zero but the energies can take arbitrary values fulfilling 
the relation £ + gi > 0. 

B. Property of the zero energy states for 2tt rotation 

In general eigenfunctions must have some definite 
properties with respect to the 2tt rotation of the az- 
imuthal angle (j) in the original (x, y) plane corresponding 
to statistical properties of the eigenstates, for instance, 
ij{9,cj) + 2n) = ^lj{9,(f>) for integer-spin states (bosonic 
states) and 4'^{9, (f> + 2n) = —'ip^{9, (f>) for ^-spin states 
(fermionic states). We see that in the latter condition 
(fermionic condition) the eigenstates first return to the 
original states by 47r rotation. It is known that states 
having new statistical properties which are different from 
Bose-Einstein and Fermi statistics appear in two dimen- 
sions and they are called anyon [45-49] . We may consider 
the eigenstates which first return to the original states by 
21^ rotation (I = integers > 3), of which condition in two 
dimensions is expressed by -0'"'((/) -I- 2tt) — e^'^^'^/^4'^{(j)) 
{k — positive integers < I) corresponding to anyon states. 
The condition may be called anyonic condition. Here we 
shall study the construction of the eigenfunctions under 
some different conditions. In this argument we represent 
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the eigenfunctions in terms of ■(/'(|(ua(Q;)) and (f)'^ {ua{ct)) 
given in (|4|) and (p^). 

(i) Cases for a = ±n with n = 1, 2, 3, • • • 

Only the normal condition for the bosonic states 
(bosonic condition) 

vl'g(p,(^(a)+2^) = *g(p,(p(a)) (18) 

are available, because (p{a) = tp + a/a and then Ua{a) do 
not change in the addition of 2a7r to their phases for the 
choices 

a = ±n with n = 1,2,3,- •• . (19) 

All the eigenfunctions ■!/'o^(Ma(Q!)) and (/),^(ua (a)) includ- 
ing the freedom of the arbitrary angle a are satisfied by 
the bosonic condition. 

(ii) Cases for a = ±{2n — l)/2 with n = 1, 2, 3, • • • 
Two types of the conditions 

%{p, ^{a) + 27r; ±) = ±%{p, ^{a); ±) (20) 

are available. Since the addition of 2a7r to the phase of 
Ua{a) changes only the signs of Ua{a) for the choices 

a = ±(2n - 1) /2 with n = 1, 2, 3, • • • , 

we see that the linear combinations 

V'g(Ma(«);±) = {^^{Ua{a))±^4iUa{a)))/V2 

0g(w„(a);±) = (0±(wa(a))±0^rK(«)))/V2 (21) 
fulfill the ± conditions, respectively. 

(iii) Cases for rational numbers 

Let us here study the case for rational numbers writ- 
ten by irreducible fractional numbers — m/l, which 
satisfies the relations < < 1. We can introduce 
generalized conditions such that 

^^^''{p, (^(a) + 2tt; e'^fc-/') ^ ^^2kn/l^a. (^^^ (^(a); e*2fc./'), 

(22) 

where k is zero or a positive integer being less than I. 
It is transparent that the cases for (I — m — 1) and 
{I = 2,m = 1), respectively, correspond to the cases of 
(i) and (ii). The eigenfunctions satisfying the conditions 
are obtained as 

/-I 

V'o''K.(");%) = ^{Vkyi'oiMa + '^riarTT)), (23) 

n=0 

where < k < I and 



Since the phase rj^ does not depend on the numerator 
m of the rational number = m/l, we see that the 
eigenfunctions for all with the same denominator I can 
be expressed by the functions with the same property for 
the rotation of the angle 27r in the {x, y) plane. We can, of 
course, describe the eigenfunctions for the case with ga < 
by replacing ip^ {ua{ct+2narTr)) with (j)^ {ua{a+2narn)) 
in the right-hand side of (p3|). 

Since all the rational numbers q are expressed by g = 
Ur + n with n — 0,±1,±2, • • • and the addition of the 
integers to does not change the above argument at all, 
the eigenfunctions for q satisfying the anyonic conditions 
are obtained by replacing with q. 

At this moment it is hard to answer the question 
whether the choice of irrational numbers for a is phys- 
ically meaningful or not. 

Note here that all the solutions described by the plane 
waves fulfill the same condition with respect to the 27r 
rotation in the {ua,Va) plane, i.e., the bosonic condi- 
tion, because the addition of 27r to the angle ipa does 
not change Ua{a)- We see that in order to represent the 
whole {ua, Va) plane the double sheets of the {x, y) planes 
(Riemann surface) are needed for the choice of a = ±1/2. 
In general for the choice of = ±1/^ the I sheets of the 
(x, y) plane like I spiral sheets are required to cover the 
whole {ua,Va) plane. We may consider that the case for 
a = can be understood in the limit of Z — > oo, where 
the infinite spiral sheets are needed in the {x, y) plane. 

C. Infinite degeneracy of the zero energy states 

A kind of degeneracy arising from the angle of the 
incoming particle with respect to the x axis has been 
pointed out in section 2. We, however, see that the 
zero energy states have another type of infinite degener- 
acy that has been already solved in the two-dimensional 
PPB In the PPB the degeneracy arises from the 
pairing property of the energy eigenvalues given by 
^i(n -I- l/2)/i7, that is, the energy eigenvalues of the 
type ^i{nx — ny)h'y appear in the two-dimensional PPB 
and hence the infinitely degenerate zero energy states 
are derived for all the cases satisfying Ux = Uy. We 
see that the origin of the infinite degeneracy is due to 
the existence of the infinite number of resonances hav- 
ing the decay widths (n -I- l/2)hr) in the one-dimensional 
PPB and the coexistence of the resonance-formation and 
resonance-decay processes with equal probability in the 
two-dimensional PPB. The zero energy states are inter- 
preted as the stationary flows expressed by the incoming 
flows corresponding to the formation process and the out- 
going flows corresponding to the decay process, which 
will be shown in fig. 1 and fig. 2 of section 5. Let us see 
the degeneracy in (11) where the two-dimensional PPB 
is included. As an example we study the freedom for the 
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wavefunction 'ijj^{ua) given by ( |T^). By putting the wave- 
function /='=(wq; Wa)V'o'(^to) into (^J) where f^{u 

a'l ^a) is 

a polynomial function of Ua and Va, we obtain the equa- 
tion 

[Aa±2lka^]f^iUa;Va)=0. (24) 
OUa 

As solved in ref. a few examples of the functions / 
are given by 

fo{Ua;Va) = 1, 
fti^a^Va) = 4,kaVa, 

f^{Ua;Va)=Miklvl + l±AlkaUa). (25) 

In the two-dimensional PPB the functions are generally 
written by the multiple of the polynomials of degree n, 
H^{y/2k^x), such that 

ft{u2-.V2) = H^{^2x)-H^{^2y), (26) 

where x and y in the right-hand side should be considered 
as the functions of U2 and V2 ^ ■ Since the form of the 
equations (|2^) is the same for all a, the solutions can be 
written by the same polynomial functions that are given 
in ( |2^ ) for the PPB. That is to say, we can obtain the 
polynomials for arbitrary a by replacing U2 and V2 with 
Ua and Va i n ([2^) . Note that the polynomials H^{S) 
with ^ = ^Jmrfjhx are defined by the solutions for the 
eigenstates with E = =Fi(n -I- Xjl^hrj in one dimensional 
PPB of the type V(x) — — m7^x^/2 and they are written 
in terms of the Hermite polynomials Hn{^) as 

(^) = e±'""/4i?„ (eT»-/4^) . (27) 

(For details, see ref. |^,||.) It is remarkable that all the 
wavefunctions for arbitrary a can be represented by the 
same functions of the PPB in the {ua,Va) plane. For 
■00^(^0) we should take the polynomials f^{va',Ua) in 
which the variables Ua and Va are exchanged. 

Let us here briefly comment on the boundary con- 
ditions discussed in the last subsection. Considering 
the relations Ua{a + 2tt) = Ua{a) and Va{a + 27r) = 
Va{a), we can make the eigenstates satisfying the suit- 
able boundary conditions by using the eigenfunctions 
ipo„{ua{a)) = f+{ua(a);va{a)yip^{ua{a)) in stead of 
ipo{ua{ct)) in (3.7), ( pi|) and (p3|). Note also that the 
eigenfunctions V'o„(wa(a)) for n > 1 does not describe 
plane waves and hence they cannot be normalized in 
terms of S functions. We essentially have to treat them 
as the eigenfunctions of the conjugate space 5(M^)^ in 
GT that is expressed by 

5(K2) c £2(^2) c5(R2)\ 

where 5(R^) and £^(R^) are, respectively, Schwartz 
space and Lebesgue space in two dimensions. (For de- 
tails, see |l|,^,|[D-) We will see in the next section that 



this degeneracy plays an important role to investigate 
vortices. 

Note that we can obtain the polynomials for the wave- 
functions (f>^{ua) by replacing ka with —ika in the poly- 
nomials derived from ( p6| ) . We also easily see that in one 
dimension the equation corresponding to ( |2^ ) does not 
bring any new freedom to the plane-wave solutions. 

IV. HYDRODYNAMICAL CONSIDERATIONS 
OF THE ZERO ENERGY STATES 

In hydrodynamics conformal mappings are very power- 
ful tools to understand structures of currents. Actually 
the important hydrodynamical ideas such as the prop- 
erty of complex velocity potentials, circulations of cur- 
rents, strengths of vortices, strengths of sources and so 
forth do not change in the conformal mappings [38-41]. 
This fact means that we can simultaneously carry out 
the investigation of the hydrodynamical properties of the 
zero energy solutions for all the potentials Va{p) in the 
mapped spaces, i.e., in the {ua,Va) plane. Results for 
all the potentials with a 7^ can be obtained by the in- 
verse transformations of the conformal mappings. In this 
section we shall study the zero energy states from a hy- 
drodynamical viewpoint for the ga > cases, because the 
eigenstates for < represented by exponential grow- 
ing or dumping functions do not describe any oscillating 
waves, which will be briefly discussed in section 6. 

A. Currents and velocities 

Though states in GT are generally not normalizable, 
the probability currents are observable in physical pro- 
cesses such as in scattering processes. We shall, there- 
fore, study the currents and other quantities based on 
hydrodynamics. The probability density p{t, x, y) and 
the probability current j{t,x,y) of a state '4!{t,x,y) in 
non-relativistic quantum mechanics are defined by 

p(i,a;,y) = |^(^,a;,y)|^ (28) 
j{t, X, y) = Re [^{t, X, y)* {-ihV) i,{t, x, y)] /m. (29) 

They satisfy the equation of continuity 

|+V.,-=0. (30) 

Following the analogue of the hydrodynamical approach 
[38-41], the fluid can be represented by the density p and 
the fluid velocity v. They satisfy Euler's equation of 
continuity 

^ + V-ipv)^0. (31) 
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Comparing this equation with ( pO| ) , we are thus led to the 
following definition for the quantum velocity of a state 
tjj{t,x,y), 

_ jit,x,y) , . 

in which j{t, x, y) is given by (p9|). Notice that p and ji 
for the zero energy states do not depend on time t. 

Let us discuss them in the (ua, Va) plane. All quantities 
O defined in the (ua,fa) plane will be marked by the 
symbol 'hat' such as O that can easily be transformed 
into the quantities in the (a;, y) plane. It is apparent 
that in the (iia,fa) plane the currents of the plane waves 
tpQ {ua{a)) are represented by the same form for all a 

J"a = \Na\^hka cos a /m, 

jva = \Na\'^hka sin a /m. (33) 

Note here the following relations; 

Uaia) = UaCOSa + VaSina, Ua(0) = Wq, Va{0) ^ Va- 

When we represent the momentum in terms of the vector 
of the (ua, Va) plane as 

— [hka cos a, hka sin a) 

for ipQ {ua{c()) , the currents are generally written by 

J- = \Na?pJm. (34) 

Hence the velocities are given by 

v=Pa/m. (35) 

Following the argument of hydrodynamics (for details, 
see Appendix of ref. |^), we can introduce the complex 
velocity potential W as 



Wa = {pu^ - ipvJUI'rn- 



(36) 



(37) 



The velocity potential ^ and the stream function ^ can 
be introduced as same as those in hydrodynamics by 

where they satisfy the following relations in the {ua,Va) 
plane; 

. _ a* _ a*- ^ _ _ d'^ 

T; ~^ J — — — - . 

OUa OVa OVa OUa 

It is known that Cauchy-Riemann's equations are satis- 
fled by the velocity potential and the stream function. 

The velocities of the Ua and Va directions in the {x^y) 
plane are given by 

Vua = haVua ^Va = K^V^ , (38) 

where ha = a(u^ + ti2)('»-i)/2a^ Hydrodynamics tells us 
that Wa describes corner flows with the angle ir/a round 
the origin. For example, in the case of the PPB with 
a = 2 Q the plane waves in the {ua,Va) plane, ^'0^(7/2), 
are expressed in flg. ^ and flg. ||. 



1 


\ 










^""^^^"^""^^ 





FIG. 1. Corner flows for t/jq (ua) in two-dimensional PPB. 

V 



FIG. 2. Corner flows for -i/jg (ita) in two-dimensional PPB. 

Note that the states multuplied by the polynomials 
and /j^ of (^ also represent the corner flows with the 
angle tt/o. 



B. Vortices in the zero energy states 

In hydrodynamics vortices are very important objects. 
In quantum mechanics, since the velocity deflned by ( |3^ ) 
diverges at the zero points of the wavefunctions, the vor- 
tices generally appear at such nodal points of the wave- 
functions [17-20]. The situation is, however, not so sim- 
ple to determine the positions of vortices, because the 
vortices do not always appear at the points where the 
wavefunctions vanish, when the currents also vanish at 
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the same points. Since the zero energy states have an 
infinite degeneracy and also the freedom of the angle a, 
we will be able to create vortex patterns having arbitrary 
number of vortices at arbitrary positions. General study 



of quantized vortices is carried out in ref. 1 20 1 . We shall 
here discuss the vortex patterns in a few simple cases of 
the linear combinations in terms of the infinite degener- 
acy. 

Let us study the vortex structures appearing in the 
linear combinations of two eigenstates constructed from 
( p^ ) and (p5|). The following discussions are carried out 
in the {ua,Va) plane, because the singularities of the ve- 
locity does not change in the conformal mappings except 
the singularity of the mappings at origin for a T^positive 
integers. The general form of the linear combination of 
two states can be written as 



(39) 



where, since the two states are not normalized, the com- 
plex coefficients appearing in the linear combination are 
included in the two wavefunctions ■01 and ^2 ■ The abso- 
lute square of ^' is evaluated as 



|*P = |^i|2 + |V;2p + 2Re(V'i>2). 



(40) 



In general a component of the current of ^! is written as 



= -Re[^'*(^^Vi+^MV'2)], 



(41) 



where fj, = Ua or Va, and A^, and are complex func- 
tions defined by 

Let us study the nodal points of jvE'P, where the vortices 
appear. We have 



|^P = |V'l|2 + |V.2p+ 21^-111^21008 0, 



(43) 



where 9 denotes the phase between tpi and ■02- It is 
trivial that nodal points appear when the following two 
conditions are fulfilled; 

|^i| = |i/,2| and 008 6* = -1. (44) 

We put the first relation into (Esh and thus obtain 



1*1 



2|V'i|'(l 



(45) 



Taking account of the same relation \ipi\ — |V'2|, the cur- 
rent is written by 



m 



'^"\A„,\cos(l)A + \B^\cos(I)b){1 + COS0) 
{\A^\sin(j3A- |B^|sin0_B)sin6'], (46) 



where (j)A and (pB are, respectively, the phases of and 
Bfj^. The velocity is evaluated as 

Vf, = 7r~[(|^Ml C0S(/)A + \Bf,\ C0S(/)b) 

+ {\A^\smq^A - \B^\sm<j)B)-^^]. (47) 

-L ~\ COS U 

We see that the second term in the bracket [ ] diverges 

by de I'Hospital's theorem, when the second condition for 
the angle, cos 9 = —1, is fulfilled. Thus we can obtain the 
condition for the divergence of the velocity 



\A^\ sin0A - \B^\ sin0B ^ 0. 



(48) 



This equation means that the functions A^^ and must 
not be real and also the imaginary parts of A^ and B^ 
must not be equal at least for one of the components 
fi^ Ua and Va- 

Now we can summarize the conditions for the determi- 
nation of the vortex positions in the linear combinations 
of two wavefunctions ipi and tp2 as follows: 
(I) 

(II) 9 = (2^ — l)7r with ^=integers, {9 is the phase between 

-01 and -02), 

(III) |^^|sin0A - |S^|sin0_B 7^ 0, (^^ and are de- 
fined in (H)). 

Let us investigate the above conditions in a few simple 
examples. 

Example (i): It is trivial that any linear combinations 
composed of the wavefunctions with the lowest polyno- 
mial (25) have no vortex, because the condition (III) is 
not fulfilled whereas nodal points satisfying the condi- 
tions (I) and (II) appear in the linear combinations. 

Example (ii): The combination of the lowest polyno- 
mial and the second one such that 

VP = i^+{ua{a)) - Cf+iuaiO);Vam^4iMO)) 

has vortices at positions fulfilling the following conditions 
derived from (I) and (II); 

Va{0) = {~ir/MC\ka, 

9 + 9c = mr, (n = integers), (49) 

where 9c' is the phase of C and 

9 = ka[UaiO) - Uaia)] 

= ka[ua{0){l ~ cos a) — ■i;a(0) sin a]. (50) 

Let us examine the relations ( ^ ) in two cases for a = 1 
and 2, where C is taken to be a real number, i.e., 9c — 0. 

Case for a = 1: In this case we have Mi(0) = x and 
i'i(O) = y and then the relations are reduced to 



y=(-l)"/4|C|A:i, 
x(\ — cosa) — ysina ~ mr/ki. 



(51) 
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All vortices appear on the two lines y = ±l/4|C|fci paral- 
lel to the X axis and they are at the cross points of the two 
lines and the lines x = (titt + (— 1)" sina/4|C|)/fci(l — 
cos a) for a ^ 0. The positions of vortices for n = 
0, ±1, ±2, ±3 are presented in fig. ||, where a = tt is taken. 
This situation is quite similar to the vortices called par- 
allel vortex lines obtained in hydrodynamics. 
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1 

2 




n =odd 
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FIG. 3. Positions of vortices for n = 0, ±1, ±2, ±3 in a 
constant potential (a = 1), which are denoted by •. 

Case for a = 2 (PPB): Since the inverse transformation 
of the conformal mapping is described by the equations 
1*2(0) = x'^ — y'^ and ^2(0) = 2xy in PPB the relations 
are given by 

2xy^{-ir/4\C\k2, 

(a;^ — y^){l — cos a) — 2a;?; sin a = mr/k2. (52) 

Vortices appear at the cross points of x^ — y^ = {mr + 
(-l)"sina/4|C|)/fc2(l-cosa) and xy ^ (-l)'78|C|fc2. 
The positions of two vortices for 71 = and other four for 
n = ±1 are figured in fig. H, where a = tt is taken. 



























A) 









FIG. 4. Positions of vortices for n = 0, ±1 in PPB (a = 2), 
which are denoted by • for n = 0, o for n = 1 and for 
n — —1, respectively. 

The vortices appear at the symmetric positions with 
respect to the origin, which are described by the cross 
points of the two equations; 

x^ -y^ ^mr/2k2, xy /8\C\k2. (53) 

We can make so many variety of the vortex patterns by 
changing of the parameters, a and C, and the zero en- 
ergy states in terms of the polynomials (p5|). Here we 
stress that, as shown in the above discussions, the higher 
polynomial solutions with rt 7^ which are not described 
by the plane waves in the (u^, Va) plane play the essential 
roles to create vortices. 

Note here that the strength of vortex is characterized 
by the circulation F that is represented by the integral 
round a closed contour C encircling the vortex such that 

r = <j> v-ds (54) 

and it is quantized as 

r = 2Tilh/m, (55) 

where the circulation number I is an integer ]r^ , p0| , p3t . 
After simple but tedious calculations, we obtain that I = 
— 1 for the vortices with n = even and I = 1 for the 
vortices with n — odd. 

Before closing this section we point out the fact that 
we can realize almost all of vortex patterns because of 
the infinite degeneracy of the zero energy solutions. The 
study of the vortex patterns will be carried out by de- 
termining the parameter a (the type of potential) and 
by finding the best linear combination in terms of the 
infinitely degenerate zero-energy-states to describe the 
vortex patterns. 
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V. VORTICES IN THREE DIMENSIONS 

We shall study vortices in three dimensions. It is obvi- 
ous that the conformal mappings given in (Q) cannot ap- 
ply in three dimensions. Schrodinger equations in three 
dimensions are generally written as 

[- — (A + + Va{x, y, z)] 'ip{x, y,z) =^£ ip{x, y, z), 

(56) 

where A = /dx^ -f /dy^. The equations, however, 
can be reduced to two-dimensional ones in cases where 
potentials are separable into 2-1-1 dimensions and the 
zero energy solutions are applicable to such cases. We 
shall discuss in the following three cases. 

(I) Cases of free motions 

Let us consider the cases where the motion of one direc- 
tion (say z) is a free motion described by e'*^^^. In these 
cases the potentials can be written by two dimensional 
potentials such as Va (p) . In such cases we have 

h^A + Va{p)] Hx,y)e^^'^ ^{£- E,) ^{x,y)e^'''^ , 

(57) 

with Ez — h^k^/2m. The zero energy states in the two 
dimensions of {x,y) appear when the relation 

£-Ez^O (58) 

is fulfilled. It is trivial that the total energy £ must be 
positive because of E^ > 0. We see that all the vortex 
patterns discussed in the last section appear in three- 
dimensional phenomena where the motion of one direc- 
tion perpendicular to the vortex plane (the (x, y) plane) 
is a free motion. Vortices in a constant magnetic field of 
z direction will be one of these cases. 

(II) Cases of exponentials 

We can present another cases where the motions perpen- 
dicular to the vortex plane are described by exponentials 
such that ip{z) — e^'^- . The equation ( p6| ) is written as 

[-|^A + Va{p)] ^(x,y)e''^' - (£ + Ez) ^{x, y)e^'\ 

(59) 

When the relation 

£ + Ez = (60) 

is fulfilled, we have the same zero energy solutions and 
then we obtain the same vortex structures of the two 
dimensions. In these cases the total energy £ must ap- 
parently be zero or negative. 

(III) Cases for separable potentials 



Even if the potential of the z component is not constant 
potentials, the eigenvalue equations are separable into 
the {x, y) and the z- components in the cases with po- 
tentials such that 

V{x,y,z) = V2{x,y) + Vi{z). (61) 

When the eigenvalues of the z component are given by 
Ez^ the equations for the {x,y) component become same 
as those of the cases (I) and (II). 

In these three cases where the vortex plane {x, y) and 
the other axis {z) perpendicular to the vortex plane 
are completely separable and then the wavefunctions are 
written by the multiplicative forms such as 'i}){x,y)ip(z), 
all vortices are described by the axial type and the 
toroidal vortices do not appear j2^, because the posi- 
tions of the vortices in the {x, y) plane do not depend on 
z. 

Here we would like to note the construction of vortices 
by plane-wave solutions in the (x, y) plane. Let us put 
the functions 

Mx,y,z) = Ar^e*('==^+'=«^)(e*'=-^ or e'^^^) 

for k^,ky,kzeR (62) 

into (|56|), in which Va{x,y,z) = is taken. In this 
case, however, the relations for zero energy (|^ and ( |60| ) 
should not be taken but the different relations £^Ez > 
must be required. Taking h^k^/2m — £ ^ Ez, we have 
the solutions same as those for the constant potential 
V2 = —ga in (|6l|), where ga — h^k'^/2m and then 
k'^ + ky — k'^. This fact implies that the parallel vortices 
discussed in section 4.2 are producible from the plane 
wave and the degenerate states with non-zero energy in 
three dimensions. 

Real vortex phenomena [24-31] appear in three dimen- 
sional spaces. Some of the vortex phenomena will be 
understood in the cases discussed above. 

VI. SHORT NOTES ON ZERO ENERGY 
SOLUTIONS FOR Ga<0 

As noted in the section 3.1, we have the zero energy 
solutions (f)^{ua{a)) — AfaC*'^''""'-"' of (15). In general 
they are unnormalizable in the {x, y) plane. In some spe- 
cial cases, however, some of the four can be normalizable. 
For example, provided that the parameters a and a are 
taken so as to fulfill the relation 

cos(a(y9 — a) > for < f < 2tt, (63) 

4'^ {ua{oi)) can be normalizable. The relation can be ful- 
filled by the suitable choices of the parameters such that 
< a < 1/2 and -(1/2 - 2a)7r < a < 7r/2. There 
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are, of course, different choices, when we take the dif- 
ferent solutions from (f>^ (uaia)) . It is very hard to an- 
swer the question whether the choice of the solutions is 
physically meaningful or not. Such solutions, however, 
possibly have some meanings in phenomena limited in 
very special regions, provided that the solutions are used 
only in the limited regions and smoothly connected to 
other functions defined outside the regions. In fact the 
solutions are used for constructing the vortices from the 
plane-wave solutions in three-dimensional space. (See the 
argument of section 5.) 

Note also here that the solutions (j)^{ua{a)) = 
MaC^'^"""'-"^ have no current because they can be taken 
as real. The higher polynomial solutions with n > 2 of 
(25) or (26) can, however, have currents because they are 
generally complex. This means that we have a possibility 
for producing vortices from these solutions even if they 
will appear only in very limited regions. 



From this relation we obtain the eigenvalue £ as 

^^~8{2NTi)^' ^^^^ 

with — {rix + ny)/2. We can directly confirm the 
eigenvalues by solving the original equation for the solu- 
tions iplx.y) — i?(/9)e''^'^ {Iz =integers) that correspond 
to the symmetric solutions of the harmonic oscillator de- 
scribed by Ux — ny. We see that, provided that one 
of the eigenvalue problems can be solved, we can also 
obtain the eigenvalues of the other equation. It is inter- 
esting that harmonic oscillator (p^) and Coulomb type 
(p^^) are mapped each other by the conformal mapping 
and there is a relation between the energy eigenvalues of 
the two potentials in two dimensions. 

VIII. CONCLUDING REMARKS 



VII. REMARKS ON NON-ZERO ENERGY 
SOLUTIONS 

We shall briefly discuss the equation for non-zero en- 
ergy given by (0) 

As noted in section 2, this equation can be read as the 
equation for determining the strength of the coupling 
constant ga of the original potential Va — —a'^gaP^^"'^^^ 
for the given energy £. We shall, however, discuss it 
from a little different standpoint. If we can solve the 
eigenvalue problem for the potential of ~a^^£ p^a^ a)/a ^ 
we can obtain the eigenvalues of the original equation 
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Let us show one example for a — 1/2, where the original 
potential is written by 

Vi/2(p) = -^.9i/2^ for .gi/2 > 0. (64) 

For real and negative eigenvalues {£ < 0) the equation 
(^ can be understood as the two-dimensional harmonic 
oscillator with the spring constant k = 8\£\. The eigen- 
values of the two-dimensional harmonic oscillator are well 
known as 



En^My = (nx + riy + l)hhj, 



(65) 



where Ux and riy are zero or positive integers, and lo — 
2^2\£\/m. Thus we have the relation 

9l/2=En^,ny- (66) 



We have shown that all Schrodinger equations with the 
symmetric potentials of the type Va (p) in two dimensions 
can be reduced to the same equation with a constant po- 
tential for the zero energy eigenstates in terms of the con- 
formal mappings and the states with the zero energy are 
in the infinite degeneracy. The degeneracy becomes not 
only the origin of the huge variety of vortex patterns but 
it will possibly be a interesting tool to investigate com- 
plicated problems of surface physics including boundaries 
as well. And the idea can be extended to phenomena in 
three dimensions. Particularly this scheme will become 
a powerful tool to vortex phenomena. Actually a vor- 
tex lattice solution has been found in this scheme [ pO[ . 
We may expect that hydrodynamical approach in quan- 
tum mechaincs presened here will open many interesting 
aspects in physics such as the investigation of vortex pat- 
terns [24-34] . We also note that the eigenstates satisfying 
the anyon boundaries will be applicable to problems in 
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condensed matter physics such as vortex matters 
and quantum Hall effects [45-49]. 

It should be noticed that some kinds of equations in 
hydrodynamics [38-41 are obtainable from the original 
eigenvalue equation (|l) by changing parameters such as 
h and mass m. This means that the conformal map- 
pings (|^) are applicable to hydrodynamical problems in 
two dimensions and the infinite degeneracy can also take 
place. The analysis in terms of the functions obtained in 
this article will also become an interesting approach in 
many aspets of hydrodynamical problems. 

Finally we would like to note that the infinite degen- 
eracy of the zero energy solutions brings infinite vari- 
ety to many-body systems with a fixed energy, which 
possibly becomes an origin of a new entropy different 
from the Boltzmann entropy [42-44]. The new entropy 
has nothing to do with the determination of usual tem- 
peratures in thermal equilibrium but the new freedom 
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stored in the new entropy can be released in thermal non- 
equilibrium |4^ . These considerations will also bring a 
very new aspect in statistical mechanics extended from 
Hilbert spaces to Gel'fand triplets |42p3]. 
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